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ABSTRACT 

Context.  Accurate  positional  measurements  of  planets  and  satellites  are  used  to  improve  our  knowledge  of  their  dynamics  and  to  infer 
the  accuracy  of  planet  and  satellite  ephemerides. 

Aims.  In  the  framework  of  the  FP7  ESPaCE  project,  we  provide  the  positions  of  Mars,  Phobos,  and  Deimos  taken  with  the  U.S.  Naval 
Observatory  26-inch  refractor  during  the  1971  opposition  of  the  planet. 

Methods.  These  plates  were  measured  with  the  digitizer  of  the  Royal  Observatory  of  Belgium  and  reduced  through  an  optimal  process 
that  includes  image,  instrumental,  and  spherical  corrections  to  provide  the  most  accurate  data. 

Results.  We  compared  the  observed  positions  of  the  planet  Mars  and  its  satellites  with  the  theoretical  positions  from  INPOP  10  and 
DE430  planetary  ephemerides,  and  from  NOE  and  MAR097  satellite  ephemerides.  The  rms  residuals  in  RA  and  Dec.  of  one  position 
is  less  than  60  mas,  or  about  20  km  at  Mars.  This  accuracy  is  comparable  to  the  most  recent  CCD  observations.  Moreover,  it  shows 
that  astrometric  data  derived  from  photographic  plates  can  compete  with  those  of  old  spacecraft  (Mariner  9,  Viking  1  and  2). 
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1.  Introduction 

The  European  Satellites  Partnership  for  Computing 
Ephemerides  ESPaCE  project  aims  at  strengthening  the 
collaboration  and  at  developing  new  knowledge,  new  technol¬ 
ogy,  and  products  for  the  scientific  community  in  the  domains  of 
the  development  of  ephemerides  and  reference  systems  for  nat¬ 
ural  satellites  and  spacecraft.  Part  of  the  ESPaCE  program  aims 
at  providing  accurate  positional  measurements  of  planets  and 
satellites  that  will  be  used  to  compute  new  orbital  ephemerides. 
Since  we  had  demonstrated  that  a  precise  digitization  and  a  new 
astrometric  reduction  of  old  photographic  plates  could  provide 
very  accurate  positions  (Robert  et  al.  2011)  the  project  leaders 
chose  to  consider  such  observations  as  a  significant  task  for  the 
ESPaCE  program. 

At  the  end  of  the  1960s  a  satellite-observing  program  was 
initiated  at  the  United  States  Naval  Observatory  (USNO,  Wash¬ 
ington  DC,  USA)  to  obtain  high-precision  observations  of  the 
Martian  system.  The  observations  employed  a  long-focus  instru¬ 
ment  to  provide  astrometric  stability  and  the  option  of  a  filter 
mask  on  Mars  that  allowed  simultaneous  measurement  of  the 
primary  and  the  faint  satelites  (Pascu  1975)  Photographic  plates 
were  taken  during  the  1971  opposition  of  Mars  with  the  USNO 
26-inch  refractor  by  Pascu  (1979)  Some  of  these  plates  were 
measured  manually  and  the  resulting  positions  used  in  support  of 
Martian  space  missions,  as  well  as  to  study  the  dynamics  of  the 
satellites  (Lainey  et  al.  2007)  These  observations  were  limited  to 
the  positions  of  the  satellites  relative  to  Mars  rather  than  equa¬ 


*  USNO  retired 


torial  (RA,  Dec.)  positions  since  dense,  accurate  star  catalogs 
were  not  yet  available.  Thus  our  intention  was  to  use  the  meth¬ 
ods  developed  in  the  USNO  Galilean  plates  analysis  in 
(|2011|l  and  Robert  et  al.  (2011)  to  analyze  the  set  of  the 
Martian  observations.  It  is  noteworthy  that  the  1971  epoch  cor¬ 
responds  to  that  of  the  Mariner  9  mission  and  to  the  beginning 
of  the  Viking  program,  thus  permitting  a  comparison  with  their 
positional  accuracies. 


Robert 

USNO 


2.  Observations 

2.1.  The  astrophotographic  plates 

Eorty  years  ago,  the  satellite  motions  were  still  not  well-known. 
In  particular,  accelerations  in  their  longitudes  needed  to  be  inves¬ 
tigated,  and  it  also  appeared  necessary  to  obtain  new  astrometric 
positions  of  the  satellites  because  of  the  needs  for  new  space 
missions  under  preparation. 

Getting  accurate  positions  of  Phobos  and  Deimos  from 
ground-based  observations  has  always  been  difficult.  Their  ap¬ 
parent  magnitude  of  11.3  and  12.4,  respectively,  are  quite  faint 
compared  to  Mars,  about  -2.0.  Eurthermore,  their  orbits  are  very 
close  to  the  planet,  and  the  scattered  light  from  Mars  increases 
the  level  of  the  sky  background  and  affects  the  satellite  positions, 
systematically  shifting  them  toward  the  planet’s  center.  To  mini¬ 
mize  these  effects,  suitable  reduction  techniques  had  to  be  used. 

Photographic  plates  used  by  Dan  Pascu  were  Kodak  103aJ 
glass  plates.  The  observations  set  was  taken  within  five  nights 
from  August  7,  1971  to  September  3,  1971  (Pascu  1979)  Each 
plate  contains  two  exposures  shifted  in  the  RA  direction,  and 

Article  number,  page  1  of  4 


Article  published  by  ihUB  bciencesl  to  be  cited  as  tittD://dx.doi.ora/l  U.  1  Ufal  /Ut)U4-b36l /k!Ul  4243841 


Report  Documentation  Page 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

06  OCT  2014 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2014  to  00-00-2014 


4.  TITLE  AND  SUBTITLE 

Astrometric  observations  of  Phobos  and  Deimos  during  the  1971 


6.  AUTHOR(S) 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

United  States  Naval  Observatory  USNO„3458  Massachusetts  Ave 
NW„Washington„DC, 20392 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Astronomy  &  Astrophysics, October  6,  2014 

14.  ABSTRACT 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OE  PAGES 

19a.  NAME  OE 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

5 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


A&A  proofs:  manuscript  no.  manuscript 


Fig.  1.  Center  of  the  digitization  (negative)  of  the  USNO  Martian  plate 
n°03017.  Phobos  and  Deimos  are  displayed  in  the  white  circles.  Phobos 
is  closest  to  the  planet.  A  star  is  displayed  in  the  white  square  at  the 
bottom  of  the  image. 


eighty  plates  resulting  in  159  exposures  were  taken.  The  expo¬ 
sure  time  of  the  photographic  plates  are  10-60  seconds.  A  small, 
neutrally  transmitting  Nichrome  film  filter  was  deposited  on  a 
Schott  GG14  yellow  filter.  The  Nichrome  filter  was  placed  over 
the  planet  in  the  observations.  This  technique  allows  the  primary 
to  be  measured  without  affecting  its  resolution,  and  it  was  used 
in  combination  with  a  V  filter  to  fit  in  the  right  wavelengths  for  a 
refracting  telescope.  The  field  of  view  is  57  arcmin  on  the  x-axis 
and  43  arcmin  on  the  y-axis,  with  a  mean  angular  scale  of  20.851 
arcsec/mm  at  full  aperture. 

Forty-nine  plates  resulting  in  95  exposures  were  selected  and 
transmitted  to  ROB  to  be  digitized  (Robert  et  al.  201 1 ;  de  Cuyper 
et  al.  201 1)  We  were  able  to  make  accurate  measurements  for  87 
positions  of  Mars  and  Phobos  and  58  positions  of  Deimos.  This 
last  number  is  different  because  the  planetary  filter  obscured  the 
image  of  the  satellite  during  the  fourth  night  of  observations.  As 
an  example.  Figure [T] shows  the  center  of  the  digitized  (negative) 
USNO  Martian  plate  n°03017  of  August  13,  1971,  which  is  a 
typical  digitized  image.  Phobos  and  Deimos  are  displayed  there. 


ground  around  each  of  the  two  was  subtracted  with  a  3D  surface 
fitting  by  least  squares.  Even  if  the  planet  was  filtered  to  reduce 
its  brightness,  a  high  non-regular  sky  background  gradient  in  the 
satellites’  area  still  had  to  be  removed  (see  Figure^  to  determine 
the  most  precise  positions.  The  obtained  plate  positions  were  fi¬ 
nally  corrected  for  the  optical  distortion  introduced  by  the  ROB 
digitizer  optical  unit  during  the  digitization  (Robert  et  al.  201 1) 


2.3.  Astrometric  calibration 


The  observations  contain  two  to  nine  UCAC4  reference  stars 
(Zacharias  et  al.  2013),  and  the  astrometric  reductions  were  per¬ 


formed  using  suitable  four  constant  functional  models  (Robert 
|201  l|l  to  provide  equatorial  (RA,  Dec.)  astrometric  positions  of 
the  planet  and  its  satellites.  In  fact,  because  simultaneously  made 
cluster  Praesepe  plates  were  available  for  this  set,  we  demon¬ 
strated  that  it  was  possible  to  reduce  the  observations  with  a 
unique  scale  factor  p  for  the  isotropic  scale  and  a  unique  ori¬ 
entation  6.  We  also  were  able  to  estimate  the  contribution  of  in¬ 
strumental  effects.  That  is  why  mean  scale  p,  mean  orientation  ff, 
and  offsets  A^  and  A^  were  modeled.  Estimated  coma-magnitude 
terms  Cx  and  Cy,  object  magnitude  in,  average  object  magnitudes 
mo,  and  temperature  contribution  terms  ei,  e2,  and  63  were  also 
introduced  in  the  determination  of  the  tangential  (X,  Y)  coordi¬ 
nates  for  each  satellite  and  planet: 


'  X  -  p  cos  0  ■  .r  -  (p  H-  ei  sin  (62?  +  63))  sin  0  ■  y  -1-  Ax-h 
Cx-x-(m-mo) 

Y-  psin0  ■  X -H  (p -I- ei  sin(e2f  +  e3))cos  0  ■  y -i- A^h-  ^ 
Cy-y-(m-  mo) . 

Depending  on  the  nights  of  observation,  the  observed  sky 
fields  and  thus  the  local  star  catalog  densities  were  different. 
As  a  consequence,  the  number  of  UCAC4  reference  stars  could 
be  too  limited  to  allow  for  a  good  calibration  of  the  images. 
When  not  enough  UCAC4  referenced  stars  were  available,  we 
chose  to  combine  the  UCAC4  catalog  with  the  NOMAD  catalog 
(Zacharias  et  al.  2004)  in  order  to  provide  the  maximum  number 
of  (RA,  Dec.)  astrometric  positions  as  possible. 

Before  adjustment,  the  measured  coordinates  of  the  objects 
were  corrected  for  spherical  effects,  such  as  the  aberration  of 
light  and  the  light  deflection  (Kaplan  et  al.  1989)  the  phase  ef¬ 
fect  according  to  Lindegren’s  theory  (Lindegren  1977)  and  the 
total  atmospheric  refraction  (|Robert|201  l|l.  All  our  observations 
are  equatorial  (RA,  Dec.)  astrometric  positions  obtained  from 
tangential  {X,  Y)  coordinates  by  using  the  gnomonic  inverse  pro¬ 
jection,  and  determined  in  an  ICRS  geocentric  reference  frame 
to  be  easily  compared  with  the  most  recent  ephemerides. 


2.2.  Plate  measurement 

The  measured  (x,y)  coordinates  of  Mars,  its  satellites,  and  the 
reference  stars  were  extracted  from  the  mosaic  EITS  images  by 
a  suitable  process.  Eirst,  we  used  the  Source  Extractor  software 
(Bertin  &  Arnouts  1996)  to  create  a  list  of  objects  that  were  de¬ 
tected  on  the  images.  The  stars  assumed  to  be  present  in  the  field 
were  identified  from  the  chosen  catalogs.  On  the  other  hand,  we 
developed  special  IDE  software  to  extract  planet  and  satellite 
positions.  Because  the  observations  were  taken  at  the  opposition 
and  minimum  distance  from  Earth,  the  Mars  positions  were  de¬ 
termined  with  a  two-dimensional  limb  fitting  by  least  squares. 
We  used  the  method  described  in  Tajeddine  et  al.  2013  to  de¬ 
termine  limb  points  and  fit  the  shape  of  the  planet  with  a  con¬ 
figurable  ellipsoid.  The  satellite  positions  were  determined  with 
a  two-dimensional  Gaussian  profile  fitting  once  the  sky  back¬ 


3.  Positioning  resuits 

We  compared  the  observed  positions  of  Mars,  Phobos,  and 
Deimos  with  their  theoretical  computed  positions  given  by  the 
INPOPIO  planetary  ephemeris  (Eienga  et  al.  2011)  and  NOE 
satellite  ephemerides  (Lainey  et  al.  2007)  We  then  focused  on 
individual  observations  for  which  the  (O-C)s  of  the  planet  and 
the  satellites  were  independently  less  than  the  3cr  value  of  their 
rms  (O-C)  in  right  ascension  and  declination.  This  concerns  73 
positions  of  Mars,  71  positions  of  Phobos,  and  51  positions  of 
Deimos.  In  the  list  available  in  electronic  form  at  the  cdsQ  the 
corresponding  geocentric  observed  positions  refer  to  the  ICRE, 

'  Full  Table  is  only  available  in  electronic  form  at  the  CDS  via 
anonymous  ftp  to  cdsarc.u-strasbg.fr  (130.79.128.5)  or  via 
http : //cdsweb . u- strasbg . fr/cgi-bin/qcat? J /A+A/ 
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Table  1.  Extract  from  the  astrometric  positions  list  of  Mars,  Phobos, 
and  Deimos  available  in  electronic  form  at  the  CDS. 


Object 

Date  (TDB) 

RA  (deg.) 

Dec.  (deg.) 

Mars 

2441170.676697 

322.898652 

-21.835786 

Phobos 

2441170.676697 

322.908902 

-21.834846 

Deimos 

2441170.676697 

322.923583 

-21.832599 

and  the  mean  time  of  observation  is  given  in  Barycentric  Dynam¬ 
ical  Time  TDB.  Starting  from  the  left  hand  column  we  provide 
the  object  name,  the  mean  TDB  date  of  observation  in  Julian 
Days,  the  geocentric  observed  right  ascension,  and  declination 
in  degrees.  Table[2gives  an  extract  of  this  list. 

The  (O-C)  distributions  of  residuals  in  equatorial  right  as¬ 
cension  and  declination  are  provided  in  Figure  and  Table 
They  show  the  difference  of  (RA,  Dec.)  coordinates  for  individ¬ 
ual  planet  and  satellites,  hence  the  observed  positions  versus  po¬ 
sitions  calculated  from  the  INPOPIO  ephemeris.  Offsets  for  each 
night  are  very  small,  except  for  the  second  night  of  observation 
in  the  RA  coordinate.  Local  systematic  errors  of  the  reference 
star  catalogs  could  explain  this  bias.  In  fact,  the  epoch  difference 
between  these  1971  plates  and  the  central  epoch  of  the  reference 
stars  is  up  to  about  20  years.  Expected  systematic  errors  of  the 
reference  stars  at  a  1971  epoch  are  about  10-40  mas.  The  data 
shown  in  Figure]^ span  about  a  month,  that  is  to  say,  15°  along 
the  path  of  Mars.  Thus  the  sets  of  reference  stars  are  different. 
With  about  two  to  nine  reference  stars  per  field,  the  astromet¬ 
ric  reduction  may  be  affected  by  accidental  errors  of  individ¬ 
ual  reference  stars.  We  assume  that  the  exposure  timing  can  be 
ruled  out  because  the  dome  clock  was  calibrated  each  night  to 
the  USNO  masterclock,  and  the  camera  shutter  could  be  opened 
and  closed  on  integral  seconds  with  a  precision  of  0.2-0. 3  sec¬ 
onds.  The  time  of  mid-exposure  thus  has  a  precision  better  than 
0.5  seconds. 

The  average  (O-C)  values  for  the  observations  made  during 
the  1971  opposition  of  Mars  are  very  low  in  RA  and  Dec.  coor¬ 
dinates.  The  external  error  of  INPOPIO  for  Mars  is  supposed  to 
be  less  than  1  mas  (Fienga  et  al.  2011),  and  we  may  deduce  from 
our  results  that  this  external  error  is  realistic  for  both  coordinates. 
However,  the  order  of  magnitude  of  the  offsets  are  remarkably 
low.  Such  accuracy  was  only  reached  once  with  photographic 
plates  in  the  analysis  of  the  USNO  Galilean  observations  (|Robert| 
2011|l.  To  estimate  the  influence  of  the  planetary  ephemeris  on 
the  results,  we  computed  the  difference  of  observed  positions 
versus  positions  calculated  from  the  DE430  ephemeris  (Folkner 
et  al.  2014)  Results  are  quite  similar  to  those  calculated  from 
the  INPOPIO  ephemeris,  so  we  may  deduce  that  both  dynamical 
models  are  equally  well-constrained  for  this  period  and  system. 
The  external  error  of  DE430  for  Mars  is  also  supposed  to  be  less 
than  1  mas,  which  is  in  good  agreement  with  our  data.  We  finally 
compared  the  observed  positions  of  Phobos  and  Deimos  with 
their  theoretical  positions  given  by  the  most  recent  MAR097  JPL 
ephemerides  (Jacobson  2010)  The  average  (O-C)  differences  in 
RA  and  Dec.  coordinates  are  below  5  mas.  The  respective  rms 
(O-C)  differences  are  below  0.5  mas.  It  is  interesting  to  see  that 
these  observations  are  very  useful  for  evaluating  the  accuracy  of 
the  ephemerides  of  Mars  and  its  satellites  since  they  have  not 
been  used  for  the  fit  of  the  numerical  integrations  leading  to  IN¬ 
POPIO  and  DE430,  and  to  NOE  and  MAR097,  for  example. 

In  addition,  we  can  see  that  equatorial  ctq.  cos  a  dispersion 
values  are  higher  than  equatorial  ctj,  which  is  in  good  agree¬ 
ment  with  the  fact  that  the  apparent  motion  of  the  system  is 
essentially  on  right  ascension,  but  not  for  Deimos.  This  error 
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Fig.  2.  (RA,  Dec.)  (O-C)  according  to  the  NOE  and  INPOPIO 
ephemerides.  The  X  axis  shows  the  UTC  date  of  observation  and  Y 
axis  the  RA  and  Dec.  (O-C)s.  Red  squares  denote  Mars,  green  circles 
Phobos,  and  blue  triangles  Deimos. 


is  not  as  systematic  as  we  could  find  from  USNO  Jupiter  ob¬ 
servations  (Robert  et  al.  2011)  or  USNO  Saturn  observations 
(Pascu  &  Schmidt  1990)  and  no  causes  were  found  in  the  pho¬ 
tographic  observations  or  in  the  reduction  method.  Colas  (1992) 
found  similar  errors  in  his  sets,  but  the  question  remained.  The 
observations  used  to  fit  Phobos  and  Deimos  models  of  motion 
could  explain  this  error  because  the  satellite  models  were  mainly 
constrained  by  the  Mariner  9  (Duxbury  &  Callahan  1989)  and 
Viking  1,  2  (Duxbury  &  Callahan  1988)  spaceraft  data  before 
the  80s.  But  their  measurement  accuracies  were  limited  because 
of  uncertainties  in  the  spacecraft  positions.  Moreover,  the  model 
of  motion  of  Phobos  was  better  constrained  both  in  right  ascen¬ 
sion  and  declination  coordinates  afterward,  thanks  to  the  Phobos 
2  mission  (Kolyuka  et  al.  1991)  In  comparison,  the  Deimos’  ob¬ 
servations  with  Phobos  2  were  rather  poor.  This  could  explain 
why  the  effect  is  only  visible  with  Deimos  data  for  the  1971  pe¬ 
riod,  so  further  investigations  are  planned. 

The  key  point  is  that  the  rms  (O-C)  for  all  observations  of 
Mars  is  54.0  mas,  for  Phobos,  53.6  mas,  and  for  Deimos,  44.8 
mas.  These  average  rms  (O-C)  correspond  to  our  observation  ac¬ 
curacies  at  the  opposition  of  1971  and  can  be  compared  to  those 
of  previous  Earth-based  observations.  Colas  (1992)  made  CCD 
observations  during  the  1988  opposition  with  the  one-meter  tele¬ 
scope  at  Pic-du-Midi.  Because  he  did  not  use  a  suitable  filtering 
to  reduce  the  brightness  of  the  planet,  he  was  not  able  to  measure 
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Table  2.  Details  of  the  (RA,  Dec.)  (0-C)  in  mas  according  to  the  NOE/INPOPIO  ephemerides,  and  to  the  MAR097/INPOP10  ephemerides. 
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-1.8 
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-1.8 

56.9 

-5.3 

50.7 

Phobos 

6.0 

58.3 

-1.7 

48.6 

6.0 

58.4 

1.2 

48.9 

Deimos 

4.3 

41.9 

-3.1 

47.3 

-0.2 

41.3 

-8.2 

47.0 

the  center  of  Mars  so  that  all  his  observations  were  intersatellite 
positions.  Colas  found  an  overall  rms  (O-C)  of  60  mas,  while  our 
overall  intersatellite  rms  (O-C)  is  34  mas.  Earlier,  Pascu  (1975) 
made  his  first  photographic  observations  of  Mars  at  the  opposi¬ 
tion  of  1967  with  the  USNO  61-inch  reflector.  The  observation 
technique  was  similar  to  the  one  presented  in  this  paper.  Pascu 
used  astrophotographic  held  plates  to  calibrate  and  reduce  his 
observations  and  then  provided  spherical  (RA,  Dec.)  astrometric 
positions.  He  found  a  mean  error  of  80  mas  for  Phobos  and  100 
mas  for  Deimos. 

Several  factors  contributed  to  the  increase  in  precision  and 
the  most  important  was  undoubtedly  the  digitization  of  the  pho¬ 
tographic  plates.  Thus  the  background  halo,  in  which  Phobos 
was  embedded,  could  be  removed  and  model  centroids  applied  to 
the  images.  The  Nichrome  film  Alter  made  it  possible  to  measure 
the  position  of  the  planet’s  image  but  probably  made  little  contri¬ 
bution  to  the  precision  in  the  positions  of  the  satellites.  Another 
factor  was  the  accurate  and  dense  star  catalog  used,  which  made 
it  possible  to  obtain  precision  RA  and  Dec.  While  the  star  cata¬ 
log  also  contributed  to  the  calibration  of  the  plates,  that  contri¬ 
bution  was  more  important  for  the  large-scale  system  of  Jupiter 
and  Saturn  (|Robett|201  l|l. 

In  the  present  work,  we  found  an  overall  rms  (O-C)  of  about 
50  mas,  that  is  to  say  17  km.  This  could  Anally  be  compared  to 
spacecraft  accuracies  of  the  same  epoch.  In  fact,  the  1971  epoch 
corresponds  to  that  of  the  last  Mariner  9  Martian  mission  and  to 
the  beginning  of  the  Viking  program.  Mariner  9  went  into  orbit 
about  Mars  in  November  1971.  The  accuracy  of  its  data  set  has 
been  first  estimated  to  be  10  km  (Duxbury  &  Callahan  1989) 
Viking  1  and  2  went  into  orbit  about  Mars  in  June  and  August 
1976.  The  accuracy  of  their  data  set  has  been  first  estimated  to  be 
4  km  (Duxbury  &  Callahan  1989)  Nevertheless,  such  numbers 
do  not  consider  the  etTor  on  the  spacecraft  position  itself.  As  a 
result,  the  rms  (O-C)  of  these  spacecraft  measurements  range  be¬ 
tween  4.5  km  and  25  km  (Lainey  et  al.  2007)  Moreover,  recent 
observations  of  Deimos  by  Phobos  2  have  a  similar  error  of  14- 
16  km.  Since  spacecraft  observations  are  used  to  acquire  several 
images  on  the  same  orbit,  errors  on  their  position  induce  non- 
Gaussian  noise  due  to  the  varying  distance  to  the  imaged  moon 
and  the  sometimes  limited  number  of  observations.  As  a  con¬ 
sequence,  we  conclude  that  astrometric  data  from  photographic 
plates  with  a  new  treatment  can  provide  observations  of  Mars 
moons  similar  in  accuracy  to  old  space  data  but  with  a  much 
more  Gaussian  error. 

4.  Conclusion 

We  analyzed  a  full  series  of  photographic  plates  of  the  Martian 
satellites  taken  at  USNO  during  the  1971  opposition  of  Mars 
and  extracted  all  the  important  information  contained  in  the  plate 
data.  We  were  able  to  correct  for  instrumental  and  spherical  ef¬ 
fects  during  the  reduction,  thereby  decreasing  the  number  of  un¬ 
known  parameters.  Thanks  to  the  new  technologies  and  new  as¬ 
trometric  catalogs,  we  were  able  to  provide  astrometric  positions 
with  never  achieved  accuracy  for  Mars  moons  in  the  1970s  from 


the  ground.  We  obtained  equatorial  RA  and  Dec.  positions  of 
Phobos  and  Deimos,  allowing  us  to  deduce  positions  of  Mars  in¬ 
directly  through  accurate  ephemerides  of  the  satellites,  but  we 
were  also  able  to  obtain  positions  of  the  planet  directly.  We  ob¬ 
tained  residuals  between  a  few  tens  of  mas,  to  better  than  100 
mas.  This  provides  a  positioning  accuracy  that  is  very  close  to 
that  of  spacecrafts  from  the  same  epoch  but  a  more  numerous 
spread  on  a  longer  interval  of  time  with  much  more  Gaussian 
error. 

The  USNO  archive  contains  more  plates  taken  in  the  inter¬ 
val  1967-1998  and  the  present  work  has  demonstrated  that  a  new 
analysis  of  these  plates  will  provide  valuable  data  for  the  dynam¬ 
ics  of  the  Martian  system.  It  is  also  important  to  notice  that  many 
useful  observations  are  available  in  most  of  the  observatories  and 
national  archives.  Thus  we  have  started  to  establish  contacts  for 
scientific  purposes,  and  more  particularly,  to  improve  the  dynam¬ 
ics  of  the  planetary  systems  with  old  data,  in  the  framework  of 
the  FP7  ESPaCE  program  and  ahead  of  the  upcoming  NAROO 
projec0  of  the  IMCCE.  We  also  look  forward  to  the  arrival  of 
the  GAIA  reference  star  catalog;  reductions  of  old  observations 
will  yield  increased  accuracy  by  eliminating  any  errors  caused 
by  reference  star  positions. 
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^  See  the  NAROO  webpage  at  http :  //www .  imcce .  fr/hosted_ 
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